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Abstract 

This  paper  is  a follow-up  to  an  earlier  paper  by 
lUeisten  which  derived  asymptotic  expansions  of  integral  transforms 
of  functions  with  logarithmic  singularities.  That  result  dealt 
with  exponentially  decaying  kernels.  In  this  paper  the  results 
are  expanded  to  include  the  case  of  oscillatory  kernels  - e.g., 
Fourier  or  Hankel  transforms. 


1. 


IntroJiict  ion 


i 


i 

1 

I 


Wt>  shall  dovolop  tho  asvinpCotic  expansion  of  a elass  of  integrals 


of  the  form 


(1.1)  1(\)  = J \i(At)f(t)dt. 

0 

Kor  this  class  we  shall  assume  tliat  h(t)  is  an  "oscillatory"  kernel; 
that  is, 

N(m)  -r 

(1.2)  h(t)  exi)f  iuit'"}  );  5’  a t ""(log  t)'\  t -*■ 

„ _ „ mn 


Mere,  Re  r f and  M(m)  is  finite  for  eacli  m.  We  assume  that  h and  f 
m 

are  infinitely  differentiable  on  (0,  a).  Kurtliermore , f(t)  is  assumed 

CO 

to  vanish  "C  smoothly"  at  T < 1. 

Thus  tlie  integral  1(A)  is  one  which  miglit  arise  from  a more 
general  integral  by  applying  tlie  appropriate  van  der  Corput  (1948) 
"neutralizer”  to  isolate  the  critical  point  at  the  origin.  The  class  of 
Integrals  is  furtl\er  distinguished  by  the  nature  of  f(t)  near  the  origin, 
namely , 


(I.  )) 


f ( I ) V.  >; 

m=() 


Oni 

X C t "'(log  t) 
mn 

n=() 


Mere,  Re  r 1 N(m)  is  finite  fi'r  each  m and  the  8 's  are  anv  complex 
m mn  ' ' 

numbers.  Furthermore,  we  assume  tliat  the  asymptotic  expansion  of  .iny 
deriv.itive  of  1 is  obt.iineii  by  d i f feiv'tU  i.it  i ng  (1.1). 


0 


1 


r 


This  work  is  a cont  iiuiat  i on  ot  an  oarlior  paper  by  one  of  tlie 

antliors  (Bloistoin,  1977)  in  wliicli  b(t.)  w.is  instead  an  "exponential" 

kernel  - iu)  in  (1  ■ *1 ) replaced  by  a negative  real  number.  Unfortunately, 

the  methoil  ot  pia)of  of  that  paper  does  not  snifice  liere.  file  relevant 

literature  for  both  classes  ot  integrals  is  cited  in  that  earlii’r  p.tper 

and  will  not  be  repeated  here.  We  tlo  remark,  liowever,  th.it,  in  comparison 

to  the  e.irlier  literature,  the  distinguishing  feature  in  both  cl.isses 

of  integrals  is  that  the  coefficients  3 m.iv  be  something  other  than 

mn 

non-negat ive  integers. 

To  carry  out  the  .in.ilvsis  below,  we  sh.ill  further  assume  that 
h(t)  is  loc.illy  integr.'ible  on  (0,  ■)  ami 

( 1 . 7i ) h ( t ) = 0 ( t '* ) , t *•  0 , a • Ke  r , Ke  a - ,i  -I. 

o O 


r 


2 . Technique  of  InteRratlon. 

We  sliall  calculate  the  asymptotic  expansion  of  1(A)  by  the 
Mellin  transform  technique.  (See  Bleistein  and  Handelsman,  1975,  Chapters 
4-b.)  To  do  so,  we  define 


(2.1)  M[h(t);  z]  = I hi(t)dt  z = X + iy 

(2.2)  M[f(t);  1-z]  = J"t~''f(t)dt 


and  use  the  Me  1 1 in-I’arseval  Theorem  to  write 


(2.3 


1(A) 


1 

2n  i 


c+i‘“ 

A'^M  [h;z]  M [f;  1 


zj  dz 


We  shall  now  quote  results  about  tins  integral.  They  are  proven 
in  the  above  cited  references  and/or  in  Titchmarsh  (19A8)  and/or  in  the 
papers  by  Handelsman  and  Lew  listed  in  the  references. 

(i)  M [h;  z]  exists  and  is  analytic  for  l+a<x'^  1+r^ 

(ii)  M [li;  z]  may  be  analytically  continued  as  a holomorphic 

function  to  tlie  right  half  plane  1 + a < x,  however. 


(2.4) 

that  is  its 
(i  ii) 


rate  of  growth  on  vertical  lines  increases  with  x. 
M [f ; 1 - zl  Is  analytic  for  x < Re  a 

o 


(2.5) 


(iv)  For  the  Bromwicti  contour  in  (2.3) 

! + ;i  < c < Rt*  iTi 

o 


3 


The  asymptotic  expansion  of  I(>)  is  generated  by  "moving"  the 
Bromwich  contour  to  the  right,  thereby  picking  up  contributions  to  the 


asymptotic  expansion  from  the  singularities  of  the  analytic  continuation 
of  M [f ; 1 - z]  . To  allow  for  this  deformation  of  contour,  we  must 
Impose  conditions  on  f(t)  wliicii  will  Insure  sufficient  decay  of  its 
Mellln  transform,  thereby  compensating  for  the  growth,  of  ti[h;  z]  . 

We  now  state  Theorem  I which  concerns  M [f ; 1 - z]  . 

Theorem  I : Suppose  f(t)  locally  integrable  on  (0,  1)  with  an  expansion 


f(t)  % X E c t ™(log  t) 
„ „ mn  ' ^ 

m=0  n=0 


t 0 


where  a and  N(m)  is  finite  for  each  m. 
m 


i)  M [f ; 1-z]  is  analytic  for  x < Re  + 1 

ii)  The  analytic  continuation  of  M [f;  1 - zj  to  the  right 


takes  the  form 


(2.6) 


M [f;  1 - z]  = 


N(m)*  c e'"^™V(3 


or  \ ^ I,  + 1 - z) 

Re(a  - a ) < k 1 ' m 

m o 


mn  -i  1 ) 
N ^mn'*’^ 


c (a  + 1 - z) 

• vM  rnn  rn  _ ^ . . 

3 = -e 

mn 


+ M^(z). 


Here,  in  )'  , we  exclude  the  terms  with  3 negative  integer,  wliile, 

mn 

In  T/\  wo  Inchulo  i)nly  terms  with  a negative  intocor.  The  function 

mn 


(z)  IsSl^iialyt ic  for  x < Rea  + k and  the  result  Is  correct  for  any  k. 

o 


(ill)  M[f;l-z]  = 0(|y|  ^^)as|y|-*-°° 

Result  (i)  was  stated  earlier.  The  proof  of  (11)  follows  closely 

the  proof  of  (11)  In  Theorem  A of  the  earlier  paper  by  Blelsteln.  Details 

of  this  and  a proof  of  (ill)  are  given  in  the  appendix. 

From  Theorem  I,  we  see  negative  integer  powers  of  B , lead  to 

logarithmic  branch  points,  non-negative  Integers  powers  lead  to  poles, 

and  all  other  6 lead  to  algebraic  branch  points, 
mn 

The  principle  part  in  the  expansion  of  M [f;  1 - z]  about  such 
singularities  takes  the  following  form: 


Case  1 ; 


B = «,  > 0 

mn  — 


t (log  t) 


(z  - a - 1) 
m 


t + I 


Case  2: 


B = Z < 0 

mn 


t^'^dog  t)^">" 


(z  - « - 1) 

m 

(n  - d: 


i - 1 


log  (z  - a - 1) 
m 


Case  3: 


3 not  an  integer 
mn 


t (log  t)  + 1 _ B^n  + 1 


3 . Main  Result 

We  can  now  state  the  main  result  about  the  asymptotic  expansion 

of  1(A). 
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Theorem  2: 


Let  I(X)  (1.1)  be  an  absolutely  convergent  (perhaps  Improper) 


integral  with  f and  h locally  integrable  on  (0,'")  and  h satisfying  1.2 
and  l.A  with  + y > - I,  then  with  f neutr.=  ilzed  about  0 and  f « 0 
for  t 2;  1 , 1(\)  has  the  expansion 


(3.1) 


1(A)  = L 


N (m) 

T*  C J(a  ,A) 


Re  (a  -u  )<  k n=o 
m o 


mn  m mn 


N (m)  c 


+ E E' 

Re (a  -a  ) < k n=o 
m o 


mn 


7 K(a  ,«,A) 


(«-!).'  ' n 


-a  - 1 - k + e 

+ 0(\  ° ) , any  e > 0. 


Here,  for  each  clioice  of  n,  E*  indicates  those  mn's  for  which  (3  is 

mn 

not  a negative  integer,  while  E'  includes  exactly  those  mn's  for  which 

b = -8.,  a negative  integer, 
mn 

As  we  shall  see,  the  functions  J and  K are  related  to  an  orderly 
asymptotic  sequence  with  increasing  Re  a^.  Their  definition  is  fairly 
complicated  but  their  asymptotic  expansions  are  more  straight  forward. 
Recall  1(A)  is  given  by  (2.3)  repeated  here: 

c+i<x> 

I(X)  f dz 

> 2n  i J 

C-loo 

Also  recall  the  earlier  comments  on  m(1i,zJ  stated  after  (2.3). 

We  need  to  consider  1(A)  with  M (f ; 1 - z")  taking  the  form 
(2.6).  To  this  end  we  define 


6 


I 

1 


C + l“> 


(3.2) 

We  must 
Case  1 : 

a + 1 . 

(3.3) 

Case  2: 


(3.4) 


(3.5) 


sum  in 
(3.6) 


J(a, 6,  X) 


1 IT  6 


C - 1 CO 


-6-1 

z ) M (h  ; 20  d z 


look  at  two  separate  cases. 

6 ^ I,  a non-neeative  integer, 

mn 

In  this  case,  J(a,i,X)  is  given  as  a residue  of  the  Integral  at 


J(a,«,,X) 


1 (A) 


a + 1 


6 ^ i 

nut 


The  result  is 


“ C e^’^*^(log  X)^  ^ 

J(a,B,X)  -V  E -1— 

j=o  x'  r(l  + 6 - j) 


with 


_ M^^^[h;n  + i 

j ~ j! 


We  must  also  deal  with  integrands  arising  from  the  second 
. 1) . We  def ine 


c+i** 


e-1 


K(a,e,X)  = 


2ni 


X ^ (a  + I - z)log(z  - a-  l)M[ji;^  dz 


.-fa 
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The  result  is 


(3.7)  KCct.H.A)  -x.  A~“"‘  Z C (‘^■!■■^)(log  A)"*^  ^ 

j=o  J ^ 


with  Cj 


determined  by 


(3.8) 


z ^ Hi  [h ; z + u + l3 


j+d-l 

z 


For  explicit  details  on  the  above  contour  integrals,  we  refer 

the  reader  to  the  earlier  paper  by  the  second  author. 

Now,  by  referring  back  to  (3.1),  let  us  comment  on  the  nature 

of  the  asymptotic  sequence.  As  the  contour  is  "moved"  right,  the 

singularities  + 1 are  encountered  beginning  with  + 1 , continuing 

with  increasing  m.  Consider  the  possibilities  for  the  contribution  from 

+ 1 . If  is  a non-negative  Integer  for  all  n,  then  we  obtain  a 

finite  expansion  in  powers  of  log  A (3.3).  We  would  then  proceed  to 

a,  + 1.  However,  if  g is  other  than  a non-negative  integer,  for 
1 on 

some  n,  we  obtain  an  infinite  expansion  in  powers  of  log  A at  + 1 
(3.4).  In  this  case  it  makes  no  sense  to  proceed  to  + 1,  since 
it  is  already  of  lower  algebraic  order  in  A and  hence,  asymptotically 
zero  with  respect  to  the  sequence 


{ A ”*0  H log  ^)  ^on  ^ ) 


8 


If  B are  all  iion-necat  i ve  iiitecers,  we  have  the  case 
mn 


I ( \ ) '^  5: 

Re  (a  -u  ) < k 
m o 


N (m) 

T.  c 

mn 

n=o 


(- 


{\'^M[h;z] 


z=a  +1 
m 


Only  in  this  very  special  case  do  you  see  contributions  from  each 
singularity  as  the  contour  moves  right.  If  one  is  not  a non-negative 

integer,  you  get  "cauglit”  at  and  have  an  infinite  expansion  in 
integral  powers  of  log  A of  the  form  .1  or  K. 

We  shall  close  this  section  with  an  example.  Wo  consider  the 

integral 


(3. 10) 


i\t 


3/2 


I(\)  = I e Unt|  t dt 


Here  f(t)  is  a single  term  of  the  form  (1.3)  with 


, 3TTi/2 

(3.11)  a = 1,  B = 3/2,  C = e 

o oo  oo 


the  last  being  chosen  so  that 

3/2  I ,3/2 

(3.12)  C (Int)  = Int 

oo  ' 


is  real  and  positive  for  o < t < 1. 

Ovir  asymptotic  expansion  will  be  of  the  form  3.1  with  only 

terms  of  the  form  .i(u  , B , A)  since  B = V^.  We  have 

oo  oo  oo 
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(3.13) 


i(X)  'V  c r(B  + 1)  j(a  ,e  ,x) 

oo  oo  o oo 


A two-term  expansion  of  I will  be 


(3.  lA)  I(X)  ~ -ir(5/2) 


^ -KIor  -idoR 

° x2r(5/2)  ' x2r(3/2) 


with 


(3.15) 


where 


(3. 16) 


We  get 


(3.17) 


M^^^fh;  a + l] 
oo  -* 


C.  = 
J 


h(t)  = e'^^  , a = 1 
oo 


= -1  Cj  = -(1  - y)  - l(n/2) 


where  y = .57721,  the  Euler-Mascheroni  constant. 

We  have  the  following  two-term  expansion 


(3. 18) 


I(X)  'V  - 1 (.42279  + 1.570791) 

X 2 A 


1/2 


We  use  two-terms  in  order  to  see  the  imaginary  part  of  the  expansion. 
Note  from  3.14  the  error  is  oj(log  X) 

In  Table  I we  tabulate  I(X)  for  X=  10,  50,  100  and  compare 
it  with  results  for  the  real  and  imaginary  parts  of  I(X)  obtained 


10 


by  numerical  integration  using  Simpson's  Rule.  We  tabulate  log  \ as  well, 

because  it  is  the  "large"  parameter  in  the  asymptotic  expansion.  We 

-2 

also  Include  (log  A)  to  give  an  indication  of  the  percentage  error  to 

3/2 

be  expected  from  a two-term  expansion  with  leading  order  (log  ^ and 

\ / 2 

error  term  o|(log  A)  j. 
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r 
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TABLE 


A 

10 

log  A 

2.3025 

(log 

. 189 

ASYMPTOTIC  RESULT 

Real  Part 

-.0252 

Imaginary  Part 

.03575 

VIA  SIMPSON'S  RULE 

Real  Part 

-.01647 

Imaginary  Part 

.03199 

RELATIVE  ERRORS 

Real  Part 

34% 

Imaginary  Part 

10.5% 

1 


50  100 

3.9120  4.6051 

.065  .047 


-.002593  -.00085216 

•00186  .00050562 


-.002497  -.000825 

•00183  .000496 


3.7%  3.1% 

1.6%  1.9% 


Appendix 


Firstly,  let  us  comment  that  the  proof  of  (ii)  Theorem  1 follows 
almost  exactly  the  proof  of  (ii)  Theorem  4 in  the  earlier  paper  by 
Bleistein.  Following  the  notation  of  that  proof  we  now  break  Sj^(t)  into 
only  two  sums 

(1)  (2) 

(Al)  S|^(t)  = (t)  + (t) 


Here  each  of  the  function  satisfies 

(j) 

(A2)  S,  (t)  =0  for  t > 1 

k — 


and  for  t < 1 


(A3) 


(A4) 


(2) 


Re  (a  -a  ) < k 
m o 


N(m)  3 u , 

E C (log  t)  t (1  - t ) 

— - A 

n=0 


3 ^ t 

mn 


E N(m)  3a  L 

Re(a  -a  ) < k EC  (log  t)  ’’’^t  "'  ( 1 - t ) 

mo  _ mn 

n=0 

3 = - H 

mn 


Now  the  details  of  the  proof  follow.  The  proof  also  shows  the  form  of 
the  principle  parts  of  M [If ; 1 - ^ which  are  given  in  Section  2. 

Now  we  prove  (iii)  Theorem  1. 

For  M [f ; 1 - ^ we  consider  integrals  of  the  form 


3 a -z 
. , , mn  m 

(log  t)  t 


(1  - t ) dt 


o 


(A5) 


1(7.) 


r 


With  the  change  of  variable  log  t = - T,  this  becomes 


iTiB  6 + 1 - z)  ,1. 

, , , mn  mn  m . , -tk,  , 

I(z)=e  T e (1-e  )Jt 


11  we  do  (j.  + B - O integrations  by  parts,  we  have 


I(z)  = 


(I  + u - z ) 
m 


Ii+ti  -i] 

^ mr»  —* 


(Note  derivatives  of  t 


van  ish . ) 


r -tk] 

[1-e  J guarantee 


all  the  boundary  terms 


We  consider 


D>B  -0 

p -X  ( 1+U  -7.)  / \ i B 

I(z)  = e ■" 


(f)  i’ ""  [■ 


-ik  Wx 


try 

e e 


-i(l+u  -X) 


[l+B  - 0 


with  z = X + i y . 

According  to  exteiuled  R iemann-Lebesciue  Theorem  by  Olver  (p.  71)  we 
have  l(z)  = 0(11  since  the  integral  converges  uniformly  at  0 and  for 
suf 1 i e i ent 1 y 1 urge  y . 

So  ru)w 


(AlO)  M [f;  1 - z]  = 0(|y| 


- D - + e - l] 

^ mn  — * 
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I'his  p.ipcr  is  .1  fol  hiw-up  to  .in  c.irlii'r  ii.ipor  by  13 1 c i s t oi  n'wli  i eh  derived 
.isvmptolie  oxp.insions  ol  intogr.il  tr.insforms  of  functions  with  log.irithmie 
s i ng.ii  I ,ir  i t i os  . fh.it  result  do.ilt  with  I'xponont  i .il  1 v di’caving,  Uornols.  In 

this  p.ipor  t lie  results  .ire  ox|i.iudod  to  in.liuh'  the  e.iso  of  osc  ill.it  orv 
kornols  - o.g.,  rourii'r  or  ll.inkol  tr.instorms. 
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